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1. INTRODUCTION 

The following paper is an analysis of the trade-offs and requirements for a 
laser cube-corner reflector array to be placed aboard the GEOS-C satellite. 
The objective of the study is to develop an improved array design which will 
give stronger returns at lower elevation angles, and reduce pulse spread in 
the return. 


2. ORBIT 

For purpose of the analysis we will assume a 500 nautical mile (927 km) 
circular orbit. 


2.1 Range Vs Zenith Angle — Using the following figure 



h = 927 km (500 n. mi.) 
R e = 6370 km (3440 n. mi.) 


Therefore 
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R = /(7.297 x 10 3 ) 2 - (6.37 x 20 3 ) 2 sin 2 0 - 6.37 x 10 3 cos 6 
=/ 53.0 x 10 6 - 40.5 x 10 6 sin 2 9 - 6.37 x 10 3 cos 6 
= [/ 53.0 - 40.5 sin 2 6 - 6.37 cos G] x 10 3 km 
This equation is evaluated in the following table and plotted in Figure 1 . 


6 (Deg) 

R (km) 

0 

927 

10 

930 

20 


30 

1040 

40 



4 (m. 14 ) 

(Deg) 

0.738 

50 

0.748 

60 


70 

1.17 

80 


90 


R (km) 

R 4 (m. 2 ' 

1310 

2.94 

1570 

6.05 

1980 

15.5 

2620 

47.8 

3540 

88.0 


2.2 Velocity Abbe r ration Vs. Zenith Angle - We will first calculate the orbital 
velocity of the satellite 

Kmm mV 2 Kmm 

— - L. but ! = mg 

(R e + h) 2 ( R e + h ) R 2 

Therefore 


R e g _ 

(R e + h) 2 ( R e + h ) 


v T = 


/ 


'(R e + h) R 2 g 


(R e + h ) 2 




Re i 


R + h 


= /(6.37 x 10 6 ) 2 x 10 
V 7.297 X 10 6 


10 3 met er s/sec 


At zenith velocity abberration is 
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However, the velocity abberration is caused only by that component of the 
velocity normal to the line of site between the transmitter and cube-corner 
array. Denoting this component by V N we will try to calculate V N and 4> z vs 6. 
For a zenith pass, the situation is a rather simple planar problem. Referring 
to the following figure 



R s i n 9 

s i n y 

(R e +h) 


Therefore 



We will further define V* as the tangential satellite velocity component in the 
plane defined by the laser tracking station (A) satellite (B) and center of the 
earth (C). For all passes of interest (those passing at reasonably high elevation 
angles) V* = V T in magnitude (except where very close to zenith). Therefore 


3 



<£ = 


C 


= 49.8 x 10 -6 co s 


7 


= 49. 



6.37 x 10 6 
7. 297x 10 6 


sin 2 9 x 10 6 = 49.8 x 10 -6 /l - .7S sin 2 9 


or 


4>/4> z = /I - 0.7S sin 2 6 

which is evaluated in the following table and shown in Figure 2. 


£(Deg) 

4>/<P z 

6 (Deg) 

0/<£ z 

0 

1.00 

50 


10 


60 

0.66 

20 


70 


30 

0.90 

80 


40 


90 

0.5 


Since 7 is also of interest we will give the following relationship 


R "I 

e 

s i n 9 — 

’ 6.37 X 10 6 " 

R + h 

•— e 


_ 7. 29 7 x 10 6 _ 


= 0.875 sin d 


which means that V max = sin -1 0.875 = 61.0 deg. Figure 3 shows 7 vs 9 


3. SIGNAL STRENGTH VS ZENITH ANGLE FOR PERFECT ARRAY 

Now that we have the orbital parameters of importance let us calculate 
the way the signal strength of the return from the array should vary with 
zenith angle. In order to simplify the evaluation we will normalize by dividing 
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by the signal strength at zenith and denote this normalized power by W N . We 
will also assume for the present, that there is no velocity abberration, and that 
therefore the cube-corners have been made diffraction limited. The perfect 
array is approximated reasonably well by a sphere coated with cube corners to 
provide a retrodirector which has a constant specific intensity independent of 
incidence angle. Under these idealized conditions the factors which affect us 
are the range of the target (R 4 )* and the variation of atmospheric transmission 
with zenith angle. 

In Figure 4 the space loss due to the R 4 dependence of signal is presented as 
a function of 6. The ordinate is R 4 at 6 divided by the R 4 at 6 = 0 (zenith). 

In addition to space loss, a certain amount of power is lost due to scatter 
and absorption in the atmosphere. Assuming that at zenith the absorption on a 
one way pass is 0.70, the transmission will be a power function of sec 6 


T = 0.70 Sec * 


For a two way pass it is 


T 2 = [0.70 seci9 ] 2 = O.7O 2sec0 


which is plotted in Figure 5. 


(Deg) 

T 2 

0 

0.49 

10 


20 

0.47 

30 


40 

0.40 


9 (Deg) 

T 2 

50 

0.33 

60 

0.24 

70 

0.12 

80 

0.016 

90 

0 


Taking the product of the R 4 and atmospheric absorption loss we get the total 
loss as a function of zenith angle. 


*See Section 8 
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d (Deg) 

r 4 t 2/r 4 t 2 

z z 

6 (Deg) 

r 4 t 2 /r 4 t 

0 

1.00 

50 

0.167 

10 

0.97 

60 

0.059 

20 

0.797 

70 

0.0125 

30 

0.58 

80 

0.0009 

40 

0.344 

90 

0 


W w = w/w z 

_ Rz T 2 
r 4 t 2 



It should be noted that the values for T 2 are not accurate beyond 80°, overesti- 
mating the strength of absorption. 


4 S SPECIFIC INTENSITIES OF CUBE CORNER ARRAYS 

Let us define the energy per unit solid angle (intensity) returned by an in- 
dividual cube comer for a unit of incident power as the specific intensity (SI). 

For a diffraction limited cube comer SI is a function of diameter squared 

SI = K x D 2 

We will not consider wavelength dependence here, because wavelength is con- 
sidered fixed at 6943A, nor will efficiency of optical coatings etc. be considered, 
but will be assumed to be part of the constant K x . When a cube comer is used 
off axis the entrance pupil becomes a function of the incidence angle a . There- 
fore the SI must be adjusted as a function of a. Because the entrance pupil area 
tends to shrink along the axis in the plane of incidence with only a small decrease 
at right angles to this direction the area of the entrance pupil may be approxi- 
mated by the on axis area multiplied by a correction factor of [l - K 2 a] . in 
order to take care a diffraction a second [1 - K 2 a] must be used to correct for 
the increased diffraction of the smaller entrance pupil. Therefore 

SI (a) = K x D 2 [1 - K 2 a] 2 

where a is the angle of evidence with respect to the axis (in radians) and K 2 = 1 
for fused silica cube comers of the type used on previous GEOS satellites. This 
relationship has been checked experimentally and found to give reasonably accu- 
rate prediction of the SI. If the cube corners are analyzed as a sphere, the 
specific intensity will be 
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S I 


=k, tf r 

Jo 


[2 77Rsin a] [1 -K 2 a] 2 Rda 


= 2 77 R 2 Kj D 2 f sin a [1 - a] 2 d a 

Jo 


['■ . 

sin a 

4 ) 


SI = 2rrK 1 R 2 D 2 | sin a [1 - 2 a + a 2 ] d a 


= 2 77 Kj R 2 D 2 





— 


*1 

M 

* l 


sin ad a - 2 

a sin a d a + 

a 2 sin a d a 


o ^ 

o J 

0 


= 2 77 Kj R 2 D 2 [{- cos a } |’ - 2 {sin a - a cos a} | J 


+ {2 a sin a + 2 cos a - a 2 cos a } 
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= 2 77 Kj R 2 D 2 [{ 1 - 0.54} - 2 {0.840 - 0 - 0.54 + 0} 

+ { 2 x .840 - 0 + 2x 0.54 - 2 - 0.54 + 0}] 

= 2 77 K x R 2 D 2 [0.46 - 2 (.30) + (1.68 + 1.08 - 2.54)] 

= 2 77 Kj R 2 D 2 [0.08] = 77 Kj R 2 D 2 [0. 16] 

For a flat array with a circular area of radius R 

S I = 77 K x R 2 D 2 [ 1 - d 2 

If the same type of cube comer is used, the flat array is 

1/0.16 = 6.25 

times more effective on axis. However, at 

[1 - a] 2 = 0.16 
[1 - a] = 0.40 

a = 0.60 rad = 34.3 deg 

t 

the spherical array becomes better than the flat. Considerable reduction in the 
number of cube corners in a spherical array can be made if only a section of the 
sphere can be used. Since only the cube comers out to= 57.3® either side of the 
line joining the source and center of the spherical array are actually contributing 
to the return, the others can be eliminated 
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If this array is attached to the GEOS- C satellite with its axis coincident and 
parallel to the spacecraft axis, it will give a good return when at zenith, but will 
degrade when off axis. Since from Figure 3 the max y is 61® the total included 
angle would have to be 114.6 + 2 (61) = 236.6®. This would make the array ap- 
pear like a sphere for all possible parts of the orbit. However, it is not neces- 
sary to go to 8 > 80® therefore the max y required is only 59° reducing the re- 
quired part of the sphere to 232.6° a slight improvement. Since the cube corners 
where a > 1 contribute very little to the array the array can further be reduced 
with little loss in signal. 

Returning to the equation for SI for the spherical array 


S I = 2 77 Kj R 2 D 2 [{- cos a} | ® - 5 - 2 {s in a - a cos a) | 

+ (2a sin a + 2 cos a - a 2 cos a} Iq' 5 ] 

= 2 77 Kj R 2 D 2 [(1 - 0.879) - 2(0.479 - 0.439 - 0 + 0) 

+ (0.479 + 1.759 - 0.220 - 0 - 2 + 0)] 

= 2 77 Kj R 2 D 2 [0.121 -0.080 + 0.018] = 2 77 Kj R 2 D 2 [0.059] 
= 77 ^ R 2 D 2 [0.118] 
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This array has only about l/4 the area but has 


0.118 

0.160 


74% 


of the SI. Therefore at only slight reduction in performance we get a large re- 
duction in the number of cube corners. Total included angle would now be only 


[0.5 + 0.5] 57.3 + 2 [59 °] = 165.3° 


If we assume that this array is placed on the bottom of the GEOS-C and we call 
the radius of the satellite R g then 


R g = R sin 


165.3 

2 


= R sin (82.9) = 0.99 R = R 


Let us assume we make the array a hemisphere 



Assume we cut out the central 60° cone from hemisphere array to allow antennas 
etc. to be placed on the earth facing side of the satellite. Then the height of the 
array would be R sin 60 = 0.867R or for a 54" diameter satellite 0.867 x 27 = 
23.4". 
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■I 


This is obviously much too large ! We will therefore now consider other types 
of arrays. 


5 0 THE TRUNCATED PYRAMID ARRAY 

From considerations of the satellite shape (An octagonal cross section) the 
possibility of placing the cube curves on eight sloped panels as shown in the 
following diagram appears feasible 



6 
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a = slope of panel with respect to satellite axis 

fi = slope of observation axis with respect to satellite axis 

y = rotation about satellite axis with 0° the plane containing the line of 
sight and the satellite axis 

This array is very similar to the earlier array used on BE- C but lacks the 
bottom panel because space must be left for other experiments. 

Using the next figure we can calculate the incidence angle as a function of 

/3 , a , and y . 



a = sin a a + cos a. a 

a x z 

sin/3 co sya x + sin/3sinya y +cos/3a z 
cos 8 - a a • ^ = sin a sin /3 cos y + cos a cos /3 


We have assigned numbers to each of the eight panels, and assume then to have 
equal numbers of cube comers. Calculations of the effective area are shown 
in the following tables, and Figures 7 and 8. 
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/3 = 77/3 sin /0 = / T /2 cos ^ - 1/2 




Panel 

cos V 

cos 0 





1 

1 

^ (/3" sin a + COS a) 





2 

/2/2 

1 //6 . \ 
— 1— Sin a + cos al 





3 

0 

I (COS a) 





4 

-/2/2 

— I - — sin a + cos a 
2 \ 2 > 

) 




5 

-1 

i (- /3* sin a + cos a) 





6 

-^2/2 

1 / /6 f 

— 1 — — - sin a + cos a 
z \ z 

) 




7 

0 

I (cos a) 





8 

/2/2 

1 (A . ^ \ 

— 1— sin a + cos a I 







cos 9 




a -* 0 

15 

30 

45 60 

75 

90 

1 

1/2 



0.965 1.00 

0.966 

0.866 

2 & 8 

1/2 



0.795 0.780 

0.720 

0.613 

3 & 7 

1/2 



0.354 0.250 

0.129 

0 

4 & 6 

1/2 



-0.28 

— 

-0.613 

5 

1/2 



— 

% 

-0.866 
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9 (Radians) 


a -* 

0 

15 

30 45 

60 

75 

90 

1 

> 1 


0.262 

0 

0.262 

0.525 

2 & 8 

> 1 


0.655 

0.672 

0.766 

0.91 

3 & 7 

> 1 


>1 

>1 

>1 

>1 

4 & 6 

> 1 


>1 



>1 

5 

> 1 


>1 



>1 



(1 

- 9 ) Radians (Neg values = 0) 




0 

15 

30 45 

60 

75 

90 

1 

0 


0.738 

1.00 

0.738 

0.475 

2 

0 


0.245 

0.328 

0.234 

0.09 

3 

0 


0 

0 

0 

0 

4 

0 


0 

0 

0 

0 

5 

0 


0 

0 

0 

0 

6 

0 


0 

0 

0 

0 

7 

0 


0 

0 

0 

0 

8 

0 


0.245 

0.328 

0.234 

0.09 

Total 

0 


1.228 

1.656 

1.206 

0.655 
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/3 = 77/4 


Panel 

7 

sin y 

cos y 

cos y 

1 


0 

0 

1 

— “ [sin a + cos a] 

2 


77/4 

/2/2 

/2/2 

/2 pf . 1 

2 l2 J 

3 


77/3 

1 

0 

/2 , 

— [cos a] ■ 

4 


377/4 

2/2 

V2/2 

i/§" r / 2 1 

— I- — sin a + cos a 

/Jr . , 

— L- sin a- + cos aj 

A 

5 


77 

0 

-1 

6 


577/4 

V2/2 

-S2/2 

— [- sin a + COS a] 
z 

7 


377/2 

-1 

0 

[cos a] 

A 

8 


777/4 

-/J/2 

/2/2 
cos 0 

/ 2 T /2 . ' *1 

— — sin a + cos a 

a — 

0 

15 

30 

45 

60 75 90 

1 

/~2/2 

0.866 

(/6 + v/2)/4 

1 

(/6 + /2)/4 0.866 /J/2 

2 & 8 

/2/2 

0.815 


(2 + 1/2 )/4 

^/3+/2)/4 0.668 1/2 

3 & 7 

/2/2 

0.680 


1/2 

/Z/4 0 

4 & 6 
5 

/2/2 

/2/2 

0.554 

• 

(2 - /Z)/4 
0 

1 1 
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radians 


a — 

0 

15 

30 

45 

60 

75 

90 

i 

0.785 

0.525 

0.28 

0 

0.28 

0.525 

0.785 

2 & 8 

0.785 

0.620 

0.53 

0.55 

0.67 

0.840 

1.05 

3 & 7 

0.785 

0.870 

0.91 

1.05 

1.2 

>1 

1.57 

4 & 6 

0.785 

0.99 

> 1 

1.42 

> 1 

>1 

2.10 

5 

0.785 

>1 

> 1 

1.57 

> 1 

>1 

2.36 



(1 - 9 ) radians (neg. values 

= 0) 



a 

0 

15 

30 

45 

60 

75 

90 

1 

0.215 

0.475 

0.72 

1 

0.72 

0.475 

0.215 

2 

0.215 

0.380 

0.47 

0.45 

0.33 

0.16 

0 

3 

0.215 

0.180 

0.09 

0 

0 

0 

0 

4 

0.215 

0.01 

0 

0 

0 

0 

0 

5 

0.215 

0 

0 

0 

0 

0 

0 

6 

0.215 

0.01 

0 

0 

0 

0 

0 

7 

0.215 

0.18 

0.09 

0 

0 

0 

0 

8 

0.215 

0.38 

0.47 

0.45 

0.33 

0.16 

0 

Total 

1.72 

1.615 

1.84 

1.90 

1.38 

0.795 

0.215 




Total 

X COS 4 a 




a — 

0 

15 

30 

45 

60 

75 

90 

cos 4 a 

1 

0.87 

0.560 

0.25 

0.063 

0.005 

0 


1.72 

1.41 

1.03 

0.475 

0.087 

0.004 

0 




8 (1- 

a) COS 4 a 




a — 

0 

15 

30 

45 

60 

75 

90 


8 

5.15 

2.12 

0.420 

0 

0 

0 
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/? = 77/6 sin /3 = l/2 cos/3 =/3/2 


Panel 


sin 7 

cos 7 

cos 0 

1 

0 

0 

1 

(sin a+ZTcos a ) 

2 

77/4 

/ 2/2 

/ 2/2 

2 sin a + /iT cos aj 

3 

77 /2 

1 

0 

■| (/3* COS a) 

4 

377/4 

/ 2/2 

-/ 2/2 

1 / . \ 
— 1 - — sin a + /3 cos a J 

5 

77 

0 

-1 

i (-sin a + /3~ cos a) 

A 

6 

577/4 

-t/2/2 

-/2/2 

1 ( /2 . /r . \ 
— — SUl a + /3 cos a 1 

7 

377/2 

-1 

0 

i (/S’ cos a ) 

A 

8 

777/4 

-i/2/2 

Z 2/2 

~ sin a + /3 COS aj 




cos <9 


a -» 

0 

15 30 

45 

60 75 90 

1 

v/3/2 

1 

0.966 

0.866 0.707 1/2 

2 & 8 

/3/2 

0.926 

0.866 

0.739 0 0.354 

3 & 7 

v/5/2 

0.750 

0.612 

0 0 0 

4 & 6 

/T /2 

0.574 

0.364 

0 0 -0.354 

5 

/3/2 

0.500 

0.266 

0 0 - 1/2 
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6 (radians) 


a -* 

0 

15 

30 

45 

60 

75 

90 

1 

0.525 

— 

0 

0.262 

0.525 

0.785 

1.05 

2 & 8 

0.525 

— 

0.384 

0.525 

0.742 

0 

1 

3 & 7 

0.525 

— 

0.726 

0.913 

0 

0 

1.57 

4 & 6 

0.525 

— 

0.960 

>1 

0 

0 

0 

5 

0.525 

— 

>1 

>1 

0 

0 

0 


{ 1 - 9 ) radians (neg. values = 0) 



0 

15 

30 

45 

60 

75 

90 

1 

0.475 

— 

1.0 

0.738 

0.475 

0.215 

0 

2 

0.475 

— 

0.616 

0.473 

0.258 

0 

0 

3 

0.475 

- 

0.274 

0.087 

0 

0 

0 

4 

0.475 

— 

0.040 

0 

0 

0 

0 

5 

0.475 

— 

0.000 

0 

0 

0 

0 

6 

0.475 

— 

0.040 

0 

0 

0 

0 

7 

0.475 

— 

0.274 

0.087 

0 

0 

0 

8 

0.475 

— 

0.616 

0.475 

0.991 

0 

0 

Total 

3.80 


2.840 

( 1 - 9 ) 

P = 0(9 = 

1.862 

|2 

: a.) 

0.991 

0.215 

0 

a -* 

0 

15 

30 

45 

60 

75 

90 

1 

1 

0.541 

0.225 

0.046 

0 

0 

0 

2 

1 

0.541 

0.225 

0.046 

0 

0 

0 

3 

1 

0.541 

0.225 

0.046 

0 

0 

0 

4 

1 

0.541 

0.225 

0.046 

0 

0 

0 

5 

1 

0.541 

0.225 

0.046 

0 

0 

0 

6 

1 

0.541 

0.225 

0.046 

0 

0 

0 

7 

1 

0.541 

0.225 

0.046 

0 

0 

0 

8 

1 

0.541 

0.225 

0.046 

0 

0 

0 


8 

4.33 

1.80 

0.368 

0.00 

0 

0 
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75 


90 


a -* 

0 

15 

30 

45 

60 

1 

0.225 


1.0 

0.541 

0.225 

2 

0.225 


0.380 

0.225 

0.066 

3 

0.225 


0.075 

0.008 

0 

4 

0.225 


0.002 

0 

0 

5 

0.225 


0.00 

0 

0 

6 

0.225 


0.002 

0 

0 

7 

0.225 


0.075 

0.008 

0 

8 

0.225 


0.380 

0.225 

0.066 


1.80 


1.914 

1.007 

0.357 




( 1 - 

9 ) 2 


a - 

0 

15 

30 

45 

60 

1 

0.046 

0.225 

0.52 

1 

0.52 

2 

0.046 

0.144 

0.22 

0.202 

0.11 

3 

0.046 

0.032 

0.008 

0 


4 

0.046 

0 


0 


5 

0.046 

0 


0 


6 

0.046 

0 


0 


7 

0.046 

0.032 

0.008 

0 


8 

0.046 

0.144 

0.22 

0.202 

0.11 


0.370 

0.571 

0.976 

1.404 

0.74 


6. MAXIMUM ALLOWABLE SIZE FOR DIFFRACTION 
LIMITED CUBE CORNER 

If we set the diameter of the cube corner so that the first null of the Airy 
disk equals the velocity abberation at zenith, we get an upper bound for the size 
of an individual cube corner: 

a R _ < — - 49.8 x 10“ 6 

BE - D 


Then 


D < 1.22 A . 1.22 X . _ 1.22 x 10~ 6 
a BE 49 .8 X 10" 6 49.8 X 10" 6 


1.71 cm = 0.675 in . D = 5/8 inch 
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But if we build into the return of the cube comer an error equal to the 
Bradley Effect or about 49.8 /arad. by building an error into one of the angles 
of the cube corner of 49.8/5 = 10 /xrad. (2 arc sec) we can compensate for the 
velocity abberration. The effect is to split the return as shown in the follow- 
ing diagram. See also Figures 9 and 10. 



One of these two intensity peaks will fall on the receiver if the direction of 
motion is in the plane determined by the vector of satellite velocity and the line 
joining transmitter to satellite. However since the satellite is not spin stabil- 
ized we will orient the cube comer errors in the rotational positions 120° apart 
to create the toroidal ring effect shown in the next figures in plan view 




Now in order for there to be no empty spots in the toroid we will set the 
requirements that the Airy Disk radius must be 1/6 the circumference of the 
velocity abberration so that the intensity patterns overlap at about one half the 
distance to the null. Therefore 


1.22 K 
D 


77 ( 2 “be) 
6 



BE 


- 1.05 a- BE 
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Therefore 


D = 0.642" still = 5/8 inch 


It should be noted that the spots are no longer circular, but tend towards ellipses 
with 2 to 1 ratio of major to minor axis. The major axis is aligned radially. 

From a cost standpoint we would like to make the cube corners as large as 
possible. Therefore if we make all the angles of the cube comer offset from 
90° by 2 sec will also get six spots in a circle 20 sec in diameter, but the spots 
will be larger. Looking at the face of the cube corner 



Each of the pie shaped sections acts as a separate source and since each is a 
smaller fraction of the cube corner diameter the diffraction will be larger 



Approximating the sector by a circle 1/2 the diameter of the cube comer, we can 
then set the size of the cube corner to 


3 " 

D = 1 .7 1 x 2 = 3.42 cm = 1.35 inches = 1 — 

8 


Thermal distortion places a practical. limit on the size of the cube corner, 
however the distortion can be tolerated near zenith where signals tend to be 
strong. For large zenith angles the effective aperture of the cube corner is re- 
duced to such a small dimension that diffraction effects are so strong that the 
thermal distortion is completely masked. 
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PRECEDING PAGE BLANK NOT FILMED 



= 0.0908 S 3 


In order to save weight, and to allow the cube corners to be packed closer 
together, the corners are cut off at a point 1/2 the distance between the point of 
the triangle and the axis of symmetry of the prism. This does not affect the 
optical performance of the cube corner on axis and has only minor effects off 
axis. The process of cutting off the edges (called optimization) leaves the cube 
corner as shown below. 



The resulting cube corner is an equilateral hexagonal entrance pupil (base) 
with a width across flats (D) equal to 2/3 S. 
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When the corners are cut off at 2/3 for optimization each corner has a base 


B =^"ri 

c 3 


M ■ 


/T 

27 


S 2 


and a depth 


L3 J l2^Jj 


/6i /r, 


Therefore the volume of each of the removed sections is 
V =Ib t 

c 3 c c 3 L 2 7 J L 6 J 486 
The volume of the optimized cube corner is 

# 

r/s 

L 27 


V 0 = V - 3V c 


/6 

486 


S 3 


= 0.0755 S 3 


Converting to terms of D (S = 3/2 D) 


V 0 = 0.0755 



0.255D 3 


For fused silica with a density p of 2.8 x 10 3 kg/m. 3 the weight is 


or 


W Q = 0.255D 3 x 2.8 x 10 3 = 715D 3 kg. (D inm.) 
W Q = 0.715D 3 grams (D in cm. ) 
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8. CALCULATION OF EXPECTED RETURN AT ZENITH 
FOR A FLAT ARRAY 


W R = 


P P 

T T 


(e t r ) : 


A cc P cc 

” (0 R) 2 


V Ar 


VL = Return signal in photoelectrons 
P T = Power transmitted = 1 joule = 3.4 x 10 18 photons 
0 = Transmitter divergence (half angle) = 5 x 10 -4 rad. 

t a = Atmospheric transmission = 70% 

P T = Optical efficiency of transmitter = 80% 
p cc = Optical efficiency of cube corner = 80% 

9 cc = Cube corner divergence (half angle) = 5 x 10~ 5 rad. 

R = Range of Satellite = 9.27 x 10 s meters. 
t r = Transmission of receiver = 0.25 
V = Quantum efficiency of receiver 0.02 
A. = Area of receiver = 0.12 meters 2 

w 3.4 x 10 18 x 0.80 x 0.49 A cc x °- 80 x 0.25 x 0.02 

W r ' x 0 . 1 2 

3.14 (5 x 10 4 x 9.27 x 10 5 ) 2 3.14 (5 x 10“ 5 x 9.27 x 10 5 ) 2 

0.425x 10 18 A cc x 4 x 10 3 x 0.12 0.54 x 10 15 A cc x 0.12 

21 50 x 10 2 3.14 (2150) 4.60 x.l0 8 


= 0.117 x 10 7 A r „ x 0.12 = 1.4 x 10 s A 

c c c c 
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If we make the array with an area of 0.28 meters 2 . 


W R = 1.4 x 10 s x 0.2 = 28,000 photo elect rons 


9. CALCULATION OF MAXIMUM WORKING ZENITH ANGLE 

Assuming we can tolerate a signal of 40 photoelectrons we can tolerate a 
degradation of 103 from the zenith case which occurs (from Figure 6) at about 
80°, for a spherical array, (or a flat array kept pointed toward the transmitter). 
Since a spherical array is only about 16% as efficient as the flat array our al- 
lowable degradation factor is only about 160, therefore 6 from Figure 6 for the 
spherical array is about 75°. 

For the truncated pyramid array from the array graphs for /3 = 45° the 10 3 
degradation point also comes at about' 75°. 

A truncated shaped surface of revolution with a cone half angle of 45° 
would be a better surface than the pyramid because it would be rotationally 
symmetric, and is nearly as simple to construct. 


10. CONCLUSIONS 

A zenith angle of 75° can be obtained with an array with an area of 
0.28 meters 2 , if the array is shaped on the surface of a 45° cone whose axis 
is coincident to the satellite axis. Since the area of the proposed array is 
0.28 meters 2 and the diameter of the satellite is about 1.37 cm (54 in.) the 
width of the array must be approximately 


0.28 = v D W.= 4.30 W 


W = M® = 0.065 m = 2.56 
4.30 


Since the cube corners will be = 3.5 cm across flats make the array look as 
follows. (See also Figure 11.) 
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Also 5D/2 = 3.5 x 2.5 = 8.75 cm slightly larger than our area based estimate. 

In this configuration each cube corner takes up a circumferential length of 
Dtan 30° + D/cos 30 = 1.732 D = 6.14 cm therefore we should be able to get 


4 x 137 77 
6. 14 


280 


in the ring. Since the area of each is = 10.5 cm 2 this gives us a total array 
area of 


280 x 10.5 = 2940 cm 2 = 0.294 m 2 

Because of the space between cube corners, this will have to be reduced 
slightly. Assuming that 270 cube corners which will provide the desired area 
of 0.28 m 2 can be fitted in the total weight of the cube corners will be 

270 x 31.4 = 8460 gr ams = 8 . 46 kg = 18 . 6 1 bs. 
an increase in weight over the GEOS II array of 14.0 lbs. 

However, the panels should not be any heavier than before so the overall 
weight of the array has pr9bably not gone up more than 50% of the GEOS II ar- 
ray, while the efficiency compared to the GEOS n array of small cubes has 
gone up 275%. 
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Figure 3. Incidence Angle Vs. Zenith Angle 
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R 4 /R 4 — % 
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9 = ZENITH ANGLE - DEGREES 

Figure 4. Normalized Space Loss 
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Figure 5. Atmospheric Transmission 
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ELEVATION ANGLE -DEGREES 

Figure 7. Signals from GEOS-C Reflector Panels Vs. Tilt Angle 
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Figure 8. Signals from GEOS-C Reflector Panels Vs. Tilt Angle 
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Figure 9. Velocity Abberration Compensation 
by Spreading of Cube Corner Return 
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SATELLITE 




Figure 10. Velocity Abberrotion Compensation 
by Splitting of Cube Garner Return 
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Figure 11. Proposed Array Configuration 
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APPENDIX A 


NORMALIZED INTENSITIES OF A CANNONBALL TYPE 
RETROFLECTOR ARRAY ' 


This calculation determines the effectiveness of a cannonball type retrore- 
flector compared to the same number (or area) of reflectors mounted on a flat 
surface 



Normalized Intensity of Array 



2 77- R 2 (1 - 6 ) 2 sin e d 0 
4tt R 5 


1 

2 



( 1 - 6) 2 s i n 6 d 6 


if (1-26+e 2 ) sin0d0=l f sin0d0-L [ 2 6 s in Q d 6 + I f 6> 2 sin 
Jq Jq Jq Jq 


6dd 


;in 1 fa- 
- 6 sin 

2 J 


r i 


5 d 5 + A 
2 


(9 2 s i n 8 d 6 


t Sin 1 - [ - 6 cos 9 + sin 6]\ + — [- 6 2 cos 6 + 2 6 sin <9 + 2 cos 9]\ 
2 0 2 0 


Sin 1 - [-cos 1 + sin 1 - 0 ]+ — [-cos 1 + 2 sin 1 + 2 cos 1 - 2 cos 0] 
2 2 
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s * n + cos 1 - sin 1 - f° S _ + sin 1 + cos 1-1 
2 2 


sin 1 . , cos 1 , sin 13 , . 

+ 2 cos 1 . - 1 = +_cos 1 - 1 

2 2 2 2 


J N = 0.422 + 1.5 (.540) - 1 = 0.422 + 0.810 - 1 = 1.232 - 1 = 0.232 


6 - 4S° = 0.785 rod 0 2 = 0.615 


T Sin4S r a a 

T , = & cos & + sin 

U N o 


+ — [- 8 2 cos 8 + 2 8 s in 8 + 2 cos 8 ] ^ 

Oi 2 o 


-±- - .78S^ t ^-o + I ; o.615>T ^ 2>T +0 _ 

4 [_ 2 2 J 2 [_ 2 2 

— 7T* i 

= .2 St/5 - + .285— + — [- 0.307 /2*+ 0.785 /2 + /2 - 2] 

2 J 2 

= .2S/2 + .143/2 - .155 /2* + . 392< ^2* + .500 /2*- 1 


= 1.142 /2 - .298 /2 - 1 = .844/2 - 1 = 1.195 - 1 = 0. 135 


6 = 28.6 = 0.500 rod 9 2 = .250 sin 9 = .479 cos 9 = .878 


47Q 1 

: 1- [- .500 (.878) + .479 - 0] + _[- .2S (.878) + 979 + 2 (.878) - 2] 

2 2 

• 24 - [- .439 + .479] + - [- .219 + .479 + 1.755 - 2] 

2 


.29 - (.040) +-(.015) = .24 - .04 + .007 = .297 - .04 


T n= - 207 
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APPENDIX B 




REFLECTIVE AREA OF A CUBE CORNER 


This calculation is to determine the area of a cube corner as a function of 
angle of incidence of the light. Figure B-l is a sketch of the dimensions of 
the cube under analysis. This cube has been optimized for maximum on-axis 
area to weight ratio by removing the corners of the cube so that the entrance 
and exit pupils are equilateral hexagons. 

The hexagonal forward perimeter of the cube corner forms the entrance 
pupil, and its image by reflection in the cube forms the exit pupil. When observed 
from the front at an angle off axis from the line of symmetry it appears as though 
it were two hexagonal stops whose axes are coincident with the linp of symmetry 
of the prism. (See Figure No. B-2) 

To the observer it appears as shown in Figure B-3. (Assuming that the cube 
corner and observer are immersed in a fluid with an optical index equal to the 
index of the cube corner.) 

The area of the clear area is ^y inspection 

A = (w tan 30°) S + 2 = S (w tan 30° + H) 


S = (w - 2D tan I) cos I = ( w - 2 tan I ] cos I 

/2 


= w (1 - y/2 tan I) cos I 


also 


u w , , A o s w , n o w(l - v/J tan I) cos I 

H - — tan 30 = - tan 30 v 


w cos I 2 


w cos I 


= ^ tan 30° (1-/2 tan I) 


Therefore 
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Figure B-l. Geometry of Cube Corner 



ENTRANCE c EXIT 

DUAL PUPIL 


Figure B-2. Entrance-Exit Pupil Geometry 
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A = w (1 - /2* tan I) wtan 30° + ^- tan 30° (1 - /2 tan I) cos I 

= w 2 cos I tan 30° (1-/2 tan I) [ -_ n 

\2 /2 

If the cube corner is made of glass, and is not immersed in a refractive 
fluid, account of the refraction at the glass air surface must be taken. To do 
this, I is changed to I’ in all parts of the equation except cos I. Therefore 


A - vi 2 cos I tan 30° (1-/2 tan I') ( — - — tan I' 

v \ 2 2 

= l£ W 2 cos I (1 - /2 tan I') (3 - /2 tan I') 

6 


where 


At zero angle of incidence this reduces to . 

A = — w 2 tan 30° = co 2 
2 2 

This area reaches zero when 

(1 - /2 tan I') = 0 

/2 tan I' = 1 
t an I' = 

/2 

I' = 35.2° 

For mirror type cube corners this is the limit of retrodirectance. For fused 
silica cube corners (N* = 1.46) the angle is extended in 

I = sin -1 (N* sin I*) = sin* 1 (1.46 x sin 35.2°) 

= sin* 1 (0.841) = 57.3° 
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